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Entanglement of spin and orbital Kondo effect is investigated on the basis of a Kondo-type 
exchange model with twofold orbital degeneracy. By using Wilson's numerical renormalization- 
group method, we examine dynamical and thermal properties respecting the difference in time- 
reversal property of multipole operators. In the presence of particle-hole symmetry, the model 
has a new non-Fermi-liquid fixed point with a fractional entropy. The spectral intensity of the 
quadrupole susceptibility diverges in the zero-frequency limit, while the dipole susceptibility 
shows a Fermi-liquid-like behavior. This is understood by mapping to the two-channel Kondo 
model, in which the dipole moment is mapped onto the operators with the scaling dimension 
A m = 1, while the quadrupole moment onto the operators with another scaling dimension 
A e = 1/2. Even for a fairly particle-hole asymmetric case with the Fermi-liquid ground state, 
the non-Fermi-liquid behavior has significant influences in electric and thermal properties. 

KEYWORDS: orbital degeneracy, time-reversal property, multipole, numerical renormalization-group method, 
multi-channel Kondo effect, non-Fermi liquid, CezLai—^Be, Lai_ x Sr^Mn03 



§1. Introduction 

Orbital dynamics has revived interest because of new 
experimental progress in heavy-fermion systems such 
as Ce^Lai-^Bg 1-4 ) and transition metal oxides such as 
Lai-ajSr^MnOa. 5-7 ' Much attention has been attracted 
to physical origin of anomalous property in magnetic, 
elastic and transport properties in these systems. Mul- 
tipole fluctuations seem to play an essential role in real- 
izing these properties. 8 ~ n ) 

In the presence of orbital degeneracy, orbital multi- 
poles as well as spin ones constitute observables. 12-15 ) 
The simplest model to describe the combined spin and 
orbital degrees of freedom is the SU(4) symmetric spin- 
orbit model. As a consequence of the entanglement of 
spin and orbital multipoles, fluctuations with the four- 
site periodicity govern the low-energy physics in one di- 
mension. 16 ~ 20 ) On the other hand, an on-site fluctua- 
tion with conduction electrons is examined by the SU(4) 
Coqblin-Schrieffer (CS) model. The Kondo effect based 
on this model is well understood at present. 21 ' 22 ) In these 
models, however, the difference between the spin and the 
orbital multipoles is hidden in the excessively high sym- 
metry. 

In order to take into account the difference, we have 
proposed 23, 24 ) that the time-reversal property of multi- 
pole operators is essential, since the distinction between 
magnetic and electric characters should be important in 
interpreting experimental results, including dilution ef- 
fects. For instance, the dipole operator couples with 
external magnetic field, while the quadrupole operator 
interacts with distortion of the lattice. 

In the impurity problem, the CS model was ex- 
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tended toward the realistic atomic-level structure by 
Hirst. 25 ' 26 ) The model was investigated by the pertur- 
bational renormalization-group analysis. 27 ) In such con- 
text, the spin-channel Kondo model was investigated by 
several authors, 28 ~ 30 ) in which the impurity possesses 
both spin and channel degrees of freedom. The present 
authors examined more general exchange interaction in 
accordance with the point-group symmetry. 23 ) A scal- 
ing analysis was adopted to investigate the time-reversal 
difference of the multipole operators in the exchange 
model with arbitrary number of orbital degeneracies. 24 ) 
We found a new class of non-Fermi-liquid fixed point in 
the presence of the particle-hole symmetry. It was argued 
that the new fixed point exhibits anomalous behaviors 
only in electric and thermal properties. 

The purpose of this paper is to perform quantitative 
investigation of the model by a non-perturbative method. 
Namely, we use Wilson's numerical renormalization- 
group (NRG) method 31 ) to derive dynamical and ther- 
mal properties. 22 ' 32 ' 33 ) In the presence of the particle- 
hole symmetry, the system in fact has a non- Fermi-liquid 
fixed point with a fractional entropy. Near the new fixed 
point, the spectral intensity of the quadrupole suscep- 
tibility diverges in the zero- frequency limit, while the 
dipole susceptibility shows a Fermi-liquid-like behavior. 
The quantative investigation reveals that for a fairly 
particle-hole asymmetric case the non-Fermi-liquid be- 
havior has a significant influence in electric and thermal 
properties in the crossover process to the Fermi-liquid 
ground state. 

The organization of the paper is as follows. In the fol- 
lowing section, the exchange model with twofold orbital 
degeneracy is explained. Then, some preliminaries to 
the NRG calculation are given. The results of dynamical 
and thermal properties are given in §3. The nature of 
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the non-Fermi-liquid fixed point is discussed in §4. The 
final section summarizes the paper. 

§2. Model 

We begin with a Kondo-type exchange model with 
twofold orbital degeneracy. Let us write the local mul- 
tiplet as \to), where a represents a Kramers pair |, J. 
with arbitrary strong spin-orbit interaction and the crys- 
talline electric field effect. The orbital pair I = 1,2 can 
be represented by a pseudospin r z = ±1. The interac- 
tion Hamiltonian is then given by 

' ' <T f +T«T y f +(TZTf+T*T})* c -<T f 



1 



-J, 



f 



(2.1) 



where the operators with suffix c (/) denote conduction- 
(/-) electron operators at the impurity site. The 
anisotropy in r-space results from the odd character of 
t v under time-reversal operation. 23 ' Here two different 
coupling constants J m and J e have been used to charac- 
terize interactions of magnetic and electric origins. The 
SU(4) CS model can be reproduced by putting J m = J e . 
The generalization to the realistic point-group symmetry 
is straightforward. 23 ' 

In order to analyze the time-reversal property conve- 
niently, we express the local operators in alternative ba- 
sis, i.e., azimuthal component m of the fictitious spin 
j = 3/2: 

(1 T, 2 |, 2 T, 1 I) - (+3/2, +1/2, -1/2, -3/2). (2.2) 

In the new expression, one can easily see that the 
quadrupole operator connects one orbital with another, 
while the dipole and the octupole operators couple differ- 
ent spins. The former has an even property under time- 
reversal operation while the latter has an odd property. 

The interaction among operators between / and con- 
duction electrons is expressed in the scalar form of the 
spherical tensors: 24 ' 



H cx = J2 J p E ^M/MT^^/l- (2-3) 

p=l q 



It is noted that the spherical symmetry SU(2) of the 
model is adopted simply for technical convenience. Here 
T q p ^ is 4 x 4 matrix representation of the q-th component 
of the spherical tensor operator with rank p. 34 ' These 
are expressed explicitly in terms of the Clebsch-Gordan 
coefficient, 



(T (p) ) r 



(-i)f 



(2p + l) 1/2 (3/2m| M ; 3/2m'), (2.4) 



and satisfies the orthonormality, 

Tr(TW[Tj')]t) = V^ 



(2.5) 



Accordingly, the conduction and / operators are repre- 
sented in the symbolic vector notations in eq. (2.3). If 
one considers the time-reversal difference in cq. (2.3), 
i.e., J m = J\ = J3 and J e = J 2 , the Hamiltonian be- 
comes equivalent to eq. (2.1). 

To apply the NRG method to the present model, 22, 31 ~ 33 ' 



the Hamiltonian has to be cast into the semi-infinite lin- 
ear chain form. With the discretization parameter A, it 
is given by 



H N = A^- 1 '/ 2 [H k + H ex ] , 



(2.6) 



N-l 



H » = E E A-" /2 (cLc n+ i m + h.c), (2.7) 



m n=0 
3 V 



p=\ q=-p 



(2.8) 



where the "site" in the linear chain is denoted by the 
label n. The localized tensor operators for conduction 
and / electrons are defined as 

4 P) = E c L(T< p) W C cw, (2.9) 

mm' 

T^ = E \M)(T q p) )MM'(M'\. (2.10) 

MM' 

As conserved quantum numbers, the spherical sym- 
metry allows us to use the total spin in addition to the 
total charge number measured from the half-filled state. 
These quantum numbers are defined as the eigenstates 
of the following operators: 

N 

$N = C nm^mm' C nm' 
1' 

+ ]T \M)S MM ,(M'\, (2.11) 
V 



n—0 mm' 



N 



MM' 



G" = EE 



m n—0 



Cnni Gri m 



(2.12) 



where the dipole moments are related to the spherical 
tensors of rank 1, 

V5T { o ] =iS z , V5T<±{ = T-US* ± iS v ). (2.13) 

v2 

§3. Results 

In the preset work, we have used A = 3 and about 
900 states have been retained at each step of the 
renormalization-group transformation. The transforma- 
tion is given by 

C Nm C N+lm + h.C.). (3.1) 

m 

We take the unit of energy as D(l + A _1 )/2 of the mod- 
ified bandwidth due to discretization. 

3.1 Characteristic energies 

Let us first discuss the characteristic energies in the 
system. According to the scaling analysis, there are two 
characteristic energies corresponding to magnetic and 
electric responses. 24 ' In order to define the characteristic 
energies explicitly, one uses dynamical susceptibilities for 
each rank of the localized tensors, 

X P M = i r dtexpM)<[TW(t), (T,W)t]), (3.2) 
Jo 



Crossover between Fermi Liquid and non-Fermi Liquid in Orbitally Degenerate Kondo Systems 



3 



where the parenthesis denotes the average over the 
ground state. Notice that Xp 1S independent of the com- 
ponent index q due to the spherical symmetry. Since 
the condition Ji = J 3 ensures \\ = X3, Xi( w ) is termed 
in the following the magnetic susceptibility and 
X2(w) the electric susceptibility Xe( w )- 

As in the standard Kondo system, the characteristic 
energy for the magnetic response is defined by the inverse 
of the static susceptibility, i.e., T m = l/%^(0). Since the 
NRG method provides us with the spectral intensity, the 
static susceptibility can be computed by means of the 
Kramers-Kronig relation. On the other hand, the static 
quadrupolar susceptibility diverges at J e — 0. Then, we 
have to choose alternative definition for the characteristic 
energy. We use T e — aw max , where w max is the peak po- 
sition of the spectral intensity. Here we have introduced 
the coefficient a so as to make the two characteristic en- 
ergies be equal in the case of J m — J e . The numerical 
calculation gives a — 3.6952. 

The characteristic energies T m and T e is shown in Fig. 
1 as a function of J e /J m - The energies are scaled by 
the characteristic energy Tcs of the CS model. The 
dashed lines show the results by the scaling theory at 
the two-loop level. The Tk's are defined by the con- 
dition J c ff = 0.5 where the energy scale of interest be- 
comes Tk- The numerical results show good agreement 
with the scaling theory. The difference between the two 
characteristic energies is conspicuous for J e /J m < 0.3, 
while the difference for J e /J m > 1 remains small. The 
saturation of T m as J e / J m — > indicates that the mag- 
netic property is not affected by the electric one for small 
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Fig. 1. The magnetic and electric characteristic energies T m 
(square) and T e (circle) as a function of the ratio of the cou- 
plings. The values are scaled by the characteristic energy of the 
Coqblin-Schrieffer model Tcs- The dashed lines are the esti- 
mated results by the scaling analysis. 



3.2 Dynamical susceptibilities 

We show first the results of the electric response. The 
spectral intensity %" of the electric excitation is shown 
in Fig. 2. The dominant intensity, which has frequencies 
uj/T e <~ 0.3 for Je/Jm > 0.1, moves to the low- frequency 



region as J e /J m decreases. In the case of J e /J m = 0, 
the intensity remains finite at zero frequency. Thus one 
can conclude that the quadrupole response exhibits a 
non-Fermi-liquid behavior with divergent static suscep- 
tibility. 




Fig. 2. The spectral intensity of the electric excitation related 
to the quadrupole response. The zero-frequency limit of the 
intensity has a finite value in the case J e = 0. 

A plot of x" / w as a function of to/T e is shown in Fig. 3. 
Reflecting divergent behavior in the case J e = 0, the in- 
tensity in the low-frequency region is extremely enhanced 
as J e /J m decreases. However, the zero-frequency limit 
of the intensity saturates to a finite value as long as the 
electric coupling is nonzero. Although the system even- 
tually flows to the Fermi- liquid fixed point with finite J e , 
the quadrupole response cannot be characterized by an 
universal function of ui/T e . 
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Fig. 3. A plot of x"/ u as a function of ui for various Je/Jm- 
The low-frequency intensity is extremely enhanced as Je/Jm de- 
creases, reflecting divergent aspect in the case J e = 0. Note that 
both linear and logarithmic scales are used in azimuthal axis. 

The magnetic spectral intensity x m (cj) for various 
Je/Jm is shown in Fig. 4. In contrast to the electric 
response, the results for any value of J e /J m are scaled 
very well by a universal function of Lu/T m with the peak 
at uj/Tm ~ 0.3. We should emphasize that the zero fre- 
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Fig. 4. The spectral intensity of the magnetic excitation related 
to the dynamical response of the dipolc and the octupole oper- 
ators. The results for various ratio of the couplings J e /Jm are 
scaled by the function of w /T m . 



A plot of Xm/ U as a function of uo/T m is shown in Fig. 

5. The peak at uj/T m ~ 0.1 is enhanced as J e /J m de- 
creases. The zero-frequency limit satisfies the Korringa- 
Shiba relation, 35 ^ 

1™A™/™ = 1, (3-3) 

within the numerical accuracy of about 10%. It means 
that the magnetic response satisfies the Fermi-liquid re- 
lation for any ratio of the couplings J e /J m , including the 
case J e = 0. The reason for this will be explained in the 
next section. It should be stressed that for J e =0a loga- 
rithmic correction, clogw, is present as shown in the in- 
set of Fig. 5. Even though the correction is rather small 
(c ~ 10~ 2 ), this is the only trace of the non- Fermi- liquid 
fixed point in the magnetic response. The magnetic re- 
sponse can be scaled mostly by the universal function of 
u>/T m , apart from the difference in the peak height. 

3.3 T -matrix 

The T-matrix of conduction electrons can be com- 
puted 36 ' by the following correlation function, 

r-OO 

2» = -i/ dte^({j m (t),jU), (3.4) 
Jo 

where the current operator j m is given by 

= 4 E ^EE( T « P) )«- 'com'WW- (3-5) 

p q m' 

Note that the current operator is independent of the wave 
vector. 

The imaginary part of the T-matrix is shown in Fig. 

6. In the CS limit, the single peak at uj/T m ~ 0.1 is 
consistent with the previous works for the CS model 21 - 1 
and the corresponding Anderson model. 22 ) As J e /J m de- 
creases, the single peak extends down to the lower fre- 



Fig. 5. A plot of x'm/ w as a function of ui for various J e / J m . The 
logarithmic correction exists in the case of J e = (inset). 
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Fig. 6. The imaginary part of the T-matrix. The energy u) is 
scaled by the magnetic characteristic energy T m . 



quency region, and simultaneously the intensity below 
the Fermi level decreases. The fiat region straddling the 
Fermi level shrinks to extent which correlates directly 
with the crossover energy to the Fermi- liquid fixed point. 
The crossover energy can easily be estimated in results of 
the entropy and the specific heat. This suggests that the 
conventional T 2 behavior of the resistivity can be seen 
only in lower temperature. If the coupling constants are 
small enough as compared with the bandwidth, the zero- 
frequency limit satisfies the Friedel sum-rule. 

It is noted that for J e /J m < 0.01, T" relative to its 
value at the Fermi level seems antisymmetric as a func- 
tion of lu. In the limit of J e = 0, on the other hand, 
the T-matrix is a symmetric function of u due to the 
particle-hole symmetry. Thus, the structure of the T- 
matrix must change discontinuously at J e — 0. 

In the case of J e = 0, the residual resistivity asso- 
ciated with the zero-frequency limit of the T-matrix is 
larger than that in the case J e /J m ^ 0. The constant 
behavior in the T-matrix is somewhat surprising in spite 
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of the singular behavior in the quadrupole susceptibility. 
Note that in the case of the conventional two-channel 
model, 37 - 1 the self-energy of the conduction electron 38 ' 39 ) 
and the single-particle excitation of the localized elec- 
tron 40 ) behaves as M 1 / 2 due to the leading irrelevant 
operator. In the present case, however, the singular 
contribution of the quadrupole moment to the current 
opeartor vanishes because J e = in eq. (3.5). 

3.4 Entropy and specific heat 

The impurity contribution to the entropy is shown in 
Fig. 7. In the case of J e = 0, the residual entropy has the 
fractional value, In \/2, which is the same as the value in 
the two-channel Kondo model. This implies that the low- 
energy spectrum is determined by the same fixed-point 
Hamiltonian as the two-channel Kondo model. Once the 
coupling J e is switched on, the residual entropy vanishes. 
However, the entropy takes the constant value over some 
temperature range if J e /J m < 0.01. The smaller J e /J m 
gives the wider range to observe the fractional entropy. 
2.0 r 
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Fig. 8. The impurity specific heat. 
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present model (solid) : J m = 0.04, J e = 
two-channel model (dotted) : K = 0.2 
— ■ — even, — • — odd 
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Fig. 7. The temperature dependence of the impurity entropy. 
The temperature is scaled by the magnetic characteristic energy. 



Fig. 9. The flow diagram of the excited energy spectrum for the 
present model (solid line) and the two-channel Kondo model 
(dotted line). 



The impurity contribution to the specific heat is shown 
in Fig. 8. The double-peak structure in the specific heat 
reflects the release of the entropy in two steps. Two 
peaks can be separable for J e /J m < 0.01. 

§4. Non-Fermi-liquid fixed point 

Now, we turn to discuss a new class of non-Fermi- 
liquid fixed point in the presence of the particle- hole sym- 
metry, i.e., J e /J m — 0. The flow diagram of the energy 
spectrum are shown in Fig. 9, as the renormalization 
step ./V proceeds. 

The solid lines show the result of the present model, 
eq. (2.6), while the dotted lines correspond to the two- 
channel Kondo model 37 - 1 with the exchange interaction, 

£ 00' 

where the i and a specify the flavor (channel) and the 
spin degrees of freedom, respectively. S is the localized 
spin of S — 1/2. The labels referring to each line denote 
the conserved quantum numbers of the present model 



(Q,J). 

It should be emphasized that the energy spectrum of 
the present model is identical to the two-channel Kondo 
model near the fixed point. The correspondence of the 
quantum numbers are listed in Table I, where the total 
charge Q, the flavor-spin F and the spin S, have been 
used for the two-channel model. 



Table I. The correspondence of the quantum numbers between 
(a) the present model and (b) the two-channel model. 



(a) (Q,S) 


(b) (Q,F,S) 


(±1,0) 


(0, 0, 1/2) 


(0, 3/2) 


(±1,1/2,0) 


(±1,2) 


(0,1,1/2), (±2,0,1/2) 


(0,3/2), (±2,3/2) 


(±1,1/2,1) 


(±1,0), (±3,0) 


(0,0,3/2) 


(±1,1), (±1,3) 


(0, 0, 1/2), (0, 1, 1/2), (±2, 1, 1/2) 



The low-energy states in two models have the same 
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degeneracy, and the spectrum of charge excitation in the 
present model has a one-to-one correspondence to the 
spectrum of spin excitation in the two-channel model. 
Thus, the effective interaction of the present model is 
equivalent to that of the two-channel Kondo model. The 
spin degrees of freedom in the two-channel Kondo model 
is attributed to the charge ones (isospin) in the present 
model. 24 ) 

We discuss now critical behavior of the dynamical sus- 
ceptibilities based on the two-channel Kondo model. Let 
us assume that all matrix elements |(0|Tq P -'|a)| 2 decrease 
as A~ ApN . A p is the scaling dimension of a boundary 
operator in the terminology of the conformal field theory. 
Then, one obtains 



XpM ~Im^ 



l<0|Tc 



(p) 



lj — E a + iS 



A -(2A P 



(4.2) 

where we have used the fact the energy scale varies as 
A -n/2 j n ^ ne renorma ii za tion evolution. Indeed, the 
numerical calculation shows the Fermi-liquid exponent 
A m = 1.00134 for the magnetic susceptibility, while the 
non- Fermi- liquid one A e = 0.500502 for the electric one. 

Corresponding to the magnetic susceptibility in the 
present model, the matrix elements are finite only with 
the excited states a = (Q,F,S) = (0,0,1/2), (0,1,1/2), 
(±2, 1, 1/2) in the two channel model. Then, the bound- 
ary operators are characterized by the quantum numbers 
(0,0,0), (0, 1,0) and (±2, 1,0), which correspond to the 
density, the flavor-spin density and the flavor-triplet pair 
operators, respectively. These operators are classified to- 
gether by the 5*0(5) representation with its dimension- 
ality d = 10. In fact, they have the scaling dimension 
A m = 1 at the fixed point of the two-channel Kondo 
model. 41 ' Similar argument leads to the conclusion that 
the boundary operators corresponding to the electric sus- 
ceptibility belong to d — 5 with A e = 1/2, which con- 
tains the flavor-spin density operator, (0,1,0) and the 
flavor-spin singlet pair operator, (±2,0,0). 

§5. Summary 

In summary, we have investigated quantitative aspects 
of dynamical and thermal properties on spin and orbital 
Kondo effect in orbitally degenerate systems. In order to 
distinguish the spin and orbital degrees of freedoms, we 
have paid attention to the time-reversal difference of the 
multipole operators. The difference is taken into account 
by distinguishing between the coupling constant J m for 
the dipole as well as the octupole operators, and another 
constant J e for the quadrupolc operator. 

By using Wilson's NRG method, we have identified 
a new class of the non-Fcrmi-liquid fixed point in the 
presence of the particle-hole symmetry. The spectral in- 
tensity of the quadrupole susceptibility diverges in the 
zero- frequency limit as expected. By contrast, the dipole 
susceptibility shows a Fcrmi-liquid-like behavior, even 
though the fixed point has the non-Fermi-liquid excita- 
tion spectrum. This is understood by mapping onto the 
two-channel Kondo model. Namely in terms of the con- 
formal field theory, correspondence with the two-channel 
Kondo model has been made as follows: the dipole mo- 



ment in the present model is mapped onto three op- 
erators with the dimension A m = 1, the density, the 
flavor-spin density and the flavor-singlet pair operators. 
On the other hand, the quadrupole moment pretends 
as two A e = 1/2 irrelevant opeartors, the flavor-spin 
density and the flavor-spin singlet pair ones. Thus, the 
Fermi-liquid-like behavior appears only in the magnetic 
response. 

Even in realistic situations with a rather large particle- 
hole asymmetry, the anomalous behavior discussed in 
this paper should appear in electric and thermal proper- 
ties in the temperature range accessible by experiments. 
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